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Q\ \ I. INTRODUCTION 

As the concept of dynamic ensemble distribution remains elusive, a way forward to the study of nonequilibrium 
(~| systems relies partly on prototypical stochastic models The identification of scaling regimes for these latter at 
O I both large time and length scales has permitted to characterize a vast amount of nonequilibrium processes in terms 
^ of universality classes jlj. However, their asymptotic dynamics can change completely in the presence of quenched 
disorder or under space-dependent external forces This is because particle motion may well remain localized 
around strong spatial heterogeneities which prevail over stochastic fluctuations 0, [1] . 
' In this context, as well as for their value in modelling real non-equilibrium behaviour, reaction-diffusion models 
c/2 , have been much studied in recent years, with such methods as real space renormalization, yielding very detailed 
analytical results at large times In the thermodynamic limit, some universal and non- universal types of behavior 
^ _ have been reported for these systems [H and related models , yet comparatively little is known on their relaxation 
^ ■ times in finite disordered samples. In this work we focus on the dynamic exponents associated to these time scales as 
^ I well as on finite size scaling aspects which ultimately will turn out to be nonuniversal. Although these processes are 
(— ( , here schematized as one dimensional annihilating random walks in a Brownian potential 4, 5] , we expect them to be 
Q ' relevant in the description of such real cases as exciton dynamics on long disordered polymers, given the experimental 
O i success of their homogeneous counterparts ff\ . 

Typically, the microscopic dynamical rules of these models involve hard-core particles which hop randomly and 
annihilate in pairs on adjacent locations. To determine the effect of quenched disorder on dynamic exponents, here we 
^ 1 allow for both varying diffusion and reaction rates in a linear chain. So we consider the following one-step processes. 
A particle at site I (Z -I- 1) , chosen from L locations, hops with rate hi {h\) to site I + 1 {I) provided it is vacant. 
I Independently, reacting attempts take place with rate Ri whenever two particles occupy sites I and l + l. As we shall 
see in Sees. Ill and IV, the nonuniversal situation alluded to above will be associated to random orientation changes 
, in diffusion biases h — h' . This scenario is illustrated in Fig.[TJ where an evolution snapshot of these stochastic rules 
T-H ' is depicted for the case of instantaneous reactions under a binary distribution of biases. 

. Despite this simplicity, such modeling entails a variety of interesting aspects which play a crucial role in the sub- 
2? ' sequent analysis. These are: the nonequilibrium character of the process, formalisable in terms of a non-Hermitian 
quantum Hamiltonian ^8^; criticality - associated with large (domain) scales, and small energies; disorder - the dis- 
^ > tributions are the source of the remarkable non universality; particle nonconservation - bringing new features to the 
non-universal critical phenomena through the interplay between particle and hole sectors. 

This work is organized as follows. In Sec. II we present a calculation approach that exploits the quantum analogy 
' and focuses on the spectrum gap (real part) of the associated non-Hermitian Hamiltonian. For numerical convenience 
the difficulties brought in by the particle nonconservation are remedied by duplicating the degrees of freedom of the 
original process. Sec. HI outlines an analytic treatment (cf. [I,[l3|)) developed from equations of motion mappings 
combined with renormalisation and transfer matrix techniques generalised to the non-conserving non-equilibrium 
nature of the system. Among the main features and procedures (briefly presented below) are: states appearing at 
quantised values of a complex phase obtained from phase steps for each bond and sign alternations at limiting phases; 
identification and aggregation of phase steps, as generalised random walks, in general biased, in the region between the 
limiting phases; this makes the determination of phases and their distributions a first passage problem jTTI| in which 
both the character (power law or stretched exponential) of the non-universal critical dynamics and in particular the 
critical exponents are determined by the effective diffusion constant and bias, inheriting their parameter-dependences 
(non-universality). Sec. IV implements numerically the ideas developed in Sec. II to check out the theoretical results 
obtained in Sec. HI. The computation of finite size scaling properties of the gap yield the dynamic exponents z, thus 
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FIG. 1: (Color online) Evolution snapshot, corresponding to a system with 400 sites and periodic boundary conditions, with 
instantaneous reaction Ri = 1 and a binary distribution of diffusion biases hi — h'l £ {0.6, —0.2} with probability p = 0.5 
(regions in gray correspond to hi — h'l = 0.6). It can be seen that particles tend to be trapped at boundaries separating positive 
bias at the left from negative ones to the right. 



providing the fundamental scaling between length and time. In addition, numerical simulations are carried out to 
corroborate the nonuniversal findings for z. Sec. V contains a summarizing discussion along with some remarks on 
extensions of this work. 



II. FERMION DOUBLING 



The common starting point for both theoretical and numerical approaches (Sees. Ill and IV respectively), is the 
representation of the evolution operator of the stochastic process, described above, as a quantum spin-i non Hermitian 
Hamiltonian which henceforth we recast as a problem of spinless fermions via a Jordan- Wigner transformation 
( JW) ■ Furthermore, these particles become actually free so long as the constraint 

Ri=hi+h'i, (1) 

is imposed for all lattice bonds [1, ■ Specifically, it can be readily checked that in terms of fermion operators C; , C] , 
the evolution operator under these conditions reads 

L L-l 

Hc^Y. + ^'i-^^ Ci + J2{RiCi Ci+i - hi Cl^Ci ~ h[ Cj Ci+i ) +Bc, (2) 
1=1 1=1 

{Hq = h'j^). According to the parity e"-^*^ = e*" '^j '^^ of the subspace considered, here 

Be = e*"-^« ( hL Cl Cl + h'^Cl Ci - RlClCi) , (3) 

is a boundary term stemming from the JW transformation under periodic boundary conditions (PBC). 

The bilinear form imposed throughout by free fermion constraints is already very suitable for theoretical analysis 
(as given in Sec. Ill), however the pairing terms of He would pose severe size restrictions in its numerical analysis, 
demanding the consideration of matrices growing exponentially with the length of the system. In this respect, notice 
that a standard Bogoliubov transformation is not effective within our disordered non Hermitian context. Rather, 
here we follow the strategy of Ref. ^] in formally related systems and proceed to double the degrees of freedom of 
the original problem by introducing a replica Hd of the Hamiltonian He- The idea is to find a simple (i.e. disorder 
independent) unitary transformation so as to cancel out all non-conserving particle terms of the replicated system 
H^^^ = He + Hd. To this end, we put forward the following fermion operators 

fi = alCi + a2Cl + blDi + b2D} , 

gi = alCi + aiC'l + blD + biDj , (4) 
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with a's and 5's coefficients such that all pairing terms are ultimately canceled in the f — g representation. In fact, 
this can be thought of as a local Bogoliubov transformation in the real space of the double system. After some simple 
but lengthy algebra, we find that these coefficients should satisfy 
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Among the various solutions we choose, 

ai = 0.2 = 03 = — a4 = 1/2 
bl = -&2 = &3 = &4 = - V2, 

which leave us with a Hamiltonian of real parameters given by 
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As before, here i?c+D parallels the boundary terms referred to in Eq. Q, which in this language take the form 
1 
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Therefore, the number of / and g particles is conserved as all pairing terms actually cancel out, provided we restrict 
attention to subspaces having (— 1)-^<^ = (—1)-'^°. 

In passing, it is instructive to check how the exact solution of the homogeneous situation is recovered by Fourier 
transforming the // and gi operators respectively to a set of new wave fermions fq, gq. To account for the parity con- 
serving terms implicated in B^.^j^, the q-momenta must be those of the set — {±tt/L, ±3tt/L, ■ ■ ■ ,±(L — 1)tt/L} 
for the even subspace, while for the odd one the g's should belong to = {±27r/L, zLAtt/L, • • • , ±(i — 2)tt/L, 0, tt}, 
say for L even. After straightforward algebraic steps, it is readily found that the ordered Hamiltonian Ho [i.e. 
Rl — R, hi 



h, and h[ — h' \/ 1 in Eq. ([7]) ] reduces to 
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Using a simple non unitary 2x2 similarity transformation to new , fcrmion like operators [l| 
system is finally diagonalized as 
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Thus, in doubling the degrees of freedom the A^ excitations of the original problem are restored [13], but now they 
appear forked in two symmetrical branches Ag,— A_g around a Fock vacuum with "energy" RL rather than with 
vanishing value, as it would correspond to the stationary state of the original stochastic He- Hence, in general the 
doubling transformation ^ should not be regarded as a similarity one, although it preserves all the anticommuta- 
tion relations which is enough to maintain the spacing of eigenlevels, and in particular the spectrum gap. In the 
homogeneous case, this scales as 1/i^, so the usual dynamic exponent z = 2 is restituted. 

For the more interesting disordered situation, if we think of • • • ,fL,gL} as the 2L single particle states of 

the replicated Hamiltonian ^^c+n i then the single excitations of Eq. ([7]) can be obtained by diagonalizing the real non 
symmetric five-diagonal matrix (with boundary conditions) 
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where TZ = -Rj . The numerical diagonaUzation of M± [17| reveals the emergence, as before, of two symmetrical 

branches {TZ + A}, {TZ — A} of single excitations around a ground state energy TZ. Denoting by A J the lowest level 
excitation in either subspace (Re AJ — min {Re A^}), the spectrum gap g of the original evolution operator ^ 
with an initially even (odd) number of reacting-diffusing particles, is finally obtained a,s g = 2 Re Aj, with a doubly 
degenerate {g = Re Aq , with A^ being non degenerate). 

We will recapitulate these ideas in the numerical diagonalizations of Sec. IV so as to account for the theoretical 
findings of the next section. 



III. PHASE STEPS AS RANDOM WALKS 



In what follows we reconsider the non-conserving fermion Hamiltonian but without fermion doubling. One of 
the eigenvalue equations associated to that operator involves only "hole" state amplitudes a; at each site I. These 
provide an inhomogeneous term in the other equation, for the "particle" states a/'. Using two-component vectors, 
the eigenvalue equation involves 2x2 matrix coefficients which depend on the excitation "energy" A, and on bond- 
dependent biased hopping and bond-independent pair annihilation rates hi,h'i,Ri, constrained as in Eq. ([1]), to keep 
the free- fermion character. The product of the associated 4x4 transfer matrices, Il^^i provides the spectrum 
through its trace, or appropriate matrix elements (depending on boundary conditions), effectively corresponding to 
quantisation of an aggregated phase. One 2x2 part of T; is the transfer matrix ti for the self-contained subspace 
(" hole sector" ) provided by the states a; . 

From the coupled eigenvalue equations one may equivalently proceed via the corresponding 4x4 transfer matrix 
Ti, or via renormalisation using a decimation method, or by mapping ratios of successive amplitudes. Here we mainly 
employ the last approach. The transfer matrix and decimation approaches are very instructive, in clarifying respec- 
tively the role of accumulating complex phases (which provide the spectrum) and the renormalization of quenched 
probability distributions of such aggregated variables, but in what follows it is sufficient and convenient to begin with 
the transfer matrix and then move over to the mapping approach. 

For orientation and to obtain certain characteristic parameters it is useful to separate the rates (and consequently 
the transfer matrix itself) into reference uniform, Z-independent, "effective medium" parts {h,h' — R^h, indicated 
by the absence of a site label), together with the further random part involving a disorder strength 5, and (binary) 
random variables Q at each bond. Then Ti becomes Ti = T -I- A;. For any finite size L, it is possible to start from T 
and build in the disorder terms Aj completely. For that, it is convenient to work in the representation diagonalising 
T. Two (i = 1, 2 say) of the four labels on its eigenvalues relate to the "hole" sector part ti of the transfer matrix. 
The "particle" sector {i = 3,4) couples to the other sector through the inhomogeneous term referred to above. The 
eigenvalues are conveniently parametrised as A1.2 = \/h' jh exp(±i(7), where q satisfies A = 1 — 2\Jh'h cosg, with a 
similar form (with A changed sign) for A3, A4. 

The formal expansion in powers of A; gives 
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1=1 n=0 



(12) 



5 



where An and B„ come respectively from the i = 1,2, and i = 3, 4 terms in the trace, and each involves sums over n 
disorder variables. For example, 

n 

An = {isr (h'/h)^ E - E c,+™, - c,+...+.„ ll'jimq), (13) 

mi>l m2>l m„>l i=l 



where rj{mq) is a periodic function involving separations rrii of the participating disorder sites and connected with 
competing aspects of the wavelength and scales set by the local disorder and bias. Bn involves a part similar to An 
and a second more complex part generated by the inhomogeneous terms. 

It is instructive to consider the approximation to the multiple sum in Eq. (|13p coming from taking all rrii as large. 
In this case the (/-dependence in the functions ri[miq) goes into an overall exponential factor, and the product of 
these factors comes out of the multiple sum as a complex phase factor exp(±iLg) also accompanied by further q- 
independent factors. The summand is then just the product of disorder variables C... and the sum can be shown to 
be the coefficient of a;" in the generating function n^Li(l + ^Ci)- Taken together with the other factors that leads to 
J2n=o being the exponential of an aggregated complex phase with real and imaginary parts 4>i ^^^d i LKc (g), 

where 



<j,i^\n{h[/hi). (14) 

For binary random rates the sum of 0/ is like a random walk. The case of most interest is when the two possibilities 
have different signs (orientations) for the bias hi — h'i , so steps occur in both directions. The same phase step (though 
aggregated in various different ways) occurs in the parts of the corresponding approximation for X]n=o coming 
from the inhomogeneous terms; the remaining "homogeneous part" has negligible phase step. 

However, in the approximation the location of states then becomes the same as for the pure system - coming 
just from iLRe(g) since the disorder dependence of the phase is now confined to the real function (/)/. So it will 
be necessary to generalise the result obtained by taking characteristic m's as large, by retaining the registration 
between local disorder and aggregating phase which is provided by the functions rjimiq). This is most easily done 
by working with the fully-disordered transfer matrix (or directly from the original equations of motion). That gives, 
as follows, information on characteristic scales, domain structures, phase accumulations, and the region of validity of 
the approximation just discussed. For the phases from X]n=o {^^& simplest illustration) only the 2x2 subpart ti 
corresponding to the hole sector of the (fully-disordered) transfer matrix T; is needed. 

The crucial aggregating phase variable is In x/, where xi = (a;_i/a;) is the ratio of complex amplitudes ai on 
successive sites. These xi ratios satisfy 

('f )=^'(i)- (15) 
From this it follows that xi = 5/(/iJ_-^/ft,/_i), where gi is such that z; = (1 — gi)~^ satisfies 

zi+i = [{hi + K)zi + h'i_^]/[hzi + h[_^] . (16) 
From this (Mobius) map, and its first iteration yielding zij^2 i it can readily be seen that when 

z< EE [1 + 0(A)] « \zi\ « %i ^ z> , (17) 

(needing A small) the variable xi = lii "walks" according to xi+i ^ Xi + l^^(^i/^!-i) while the complex phase 
step in a; is In xi = In [g;(/iJ_-^//i;_i)] ~ In (ft,J_-^/ft,/_i) = as was found in Eq. (|14p using the approximation on 

Sn=o^"- -l^sar zi = z^ (or z<), the mapped z goes to oo (or 0, then 1, giving a vanishing amplitude) and changes 
sign. This procedure thus provides the region pT)) of validity of the approximation on ^n^gAn, and generalises it: 
the behaviour at z — z^,z^ is the proper consequence of the functions 'qiniiq), replacing the sign changes given by 
exp(±iLg) in the approximate result. 

In the corresponding discussion for the structure and phases in X)n=o (coming from "particle" states, with 
complex amplitudes O;^), linearity allows a separation into two parts: (i) without, and (ii) with, the mixing to the 
states in the "hole" subspace just treated. The discussion for the first part ((i)) is similar to that just given, and leads 
to a walk step In (oj^-^/a^) = In with c/" close to 1 if y; = {1 — cl)~^ satisfies a condition like (fT7|) . with the same 
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limits z^, to leading order in A. For small A, ipi =\n yi walks with the same step as xi = In z^;, and has the same 
behaviour at the limiting values. For part (ii), the inhomogeneous term coming from the mixing to the states in the 
other subspace leads to the occurrence again of z/ and xi = (a;_i/a/), in terms of which, for z< ^ z; <C z^ (so A is 
small and xi ^ (/ij_j^/ft-/_i)) , the amplitude ratio is /a^ = 1 + X]j=o(^'/^') where 

i+i I 

Ti = hi{a+^^ - a+) = ao n 7r' + ^0 n ^ ^ 1 ' (1^) 

0=0 3=0 

with 7i = h'l/hi, and Aq — Tq — qq/xi. 

The combination of these results from (i),(ii) gives the essence of the complicated result for the approximated 
Sn=o together with its validity criterion, and again generalises it as necessary to obtain and quantize the imaginary 
part of the phase, to enumerate the states: the limiting boundary phase conditions provide the sign reversals in 
amplitude and correspond to (reflecting and absorbing) boundary conditions on the walks described above. 

Since the walk steps of xi a-nd ^pi contain disorder we need the probability distributions P{xi) ^iid P{(fi), which 
are the same for small A. The distributions depend on the boundary conditions, and are most easily treated in a 
continuum (diffusion) picture. There, / plays the role of time, and the accumulating xi or fi plays the role of spatial 
coordinate. From the walk steps described above, within the limits set by Eq. (|17p . the distributions satisfy a standard 
biased diffusion equation having diffusion constant D and bias 6 given by 

^ - ^var In {hi/h'i) , 

b = {\n hi) - {In h[) . (19) 

The boundary conditions take the biased problem to first-passage form [ll| in which many reflections of the walk 
occur before absorption is reached, with probability exponentially small in the traversal "length". The boundary 
conditions set the limits to the excursion (of the walk), making the traversal "length" 

lnz>-lnz< =ln(Ao/A), Ao-/i, (20) 

(long, as required for the continuum picture, in the critical regime of small A). Consequently, the first passage "time" 
scale £ ^ i, is proportional to exp[(|6|/£') In (Aq/A)] , which can be written in the scaling form A~^/^ with dynamic 
exponent 

z^D/\h\. (21) 

In the opposite case of large bias, applying for nearly homogeneous systems, the effects of (Re A)^/^ dominate making 
z close to 2. 

From the exponential distribution P{x) ^ Gxp[— ( | 6 |/_D) x] for x (and similarly (p), the energy distribution is 
appreciable only at x = In ^^(A). That gives the exponential energy distribution 

P(A) cx exp[- A I 6 I h'/{D h)], for A > . (22) 

It can be seen that the dynamic exponent z is non-universal because the distribution of the random rates allows both 
signs of the walk steps and is (in general) parameter-dependent. In particular, z can diverge if the bias 6 becomes zero 
(for example, at a critical concentration pc in the binary case studied below in Sec. IV). Then a stretched exponential 
scaling form emerges as follows. For the unbiased case the first passage problem has standard cx DC relation 
between the "length" and "time" scales. That makes [ln(Ao/A)]^ cx DC, consequently 

A/Ao - exp(-c/:'') , (23) 

with (3 = 1/2 and c = CqD^/'^, with cq a numerical constant of order one. In this case (where the bias vanishes) the 
Gaussian distribution for the aggregated phase produces a gap distribution with exponential tails. Consequently, for 
& — > , the ratio of rms gap to mean gap diverges. 

IV. NUMERICAL RESULTS 

To check out the above expectations, we now investigate numerically the finite size behavior of the gap of the 
original evolution operator ^ when averaged over independent disorder realizations. To this end, we diagonalize the 
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M± matrices referred to in Eq. (|lip given in Sec. II, and whose dimensions just grow linearly with the system size. 
For concreteness and comparison with the analytical results of Sec. Ill, let us consider a binary distribution of bonds 
with parameters B12 = {^1,2 , ^1,2} 7 and probability P{B) — p5b,Bi + (1 ^ p) 63.82 ■ For a given concentration p, we 
averaged spectra typically over 6000 samples of chains with lengths 50 < L < 1000. For most p the resulting average 
gaps scale as (g) oc L^^ irrespective of parity, and eventually with non-universal dynamic exponents z = z{p), as it 
was established by Eq. (f2T|) in Sec. III. This situation is exhibited in Fig. [21 where we show a characteristic case of 
instantaneous reactions under orientation changes in diffusion biases h — h' , these latter satisfying the free ferniion 
constraint ([T]) imposed throughout. 



<g> 




100 1000 
L 



FIG. 2: Average spectrum gap (g) for a binary distribution of biases {0.6, —0.2}, and instantaneous reactions {Ri,2 = 1). The 
concentration varies from top to bottom as p = 0.5, 0.75, 0.3, 0.25. Dynamic exponents z are read off from the slopes of the 
fitting lines. For displaying purposes, data have been shifted slightly upwards with respect to p = 0.25. The inset shows the 
stretched exponential scaling expected at Pc ~ 0.2263 under free boundaries. 



On approaching the critical regime defined by the unbiased diffusion equation referred to earlier on in Sec. 
finite size effects become progressively pronounced until at a critical concentration pc, given by 



III, 



Pc = 



1 - 



In 



■ i-/ti 

, hi 



-I -1 



ln(i 



h2 



^1 > O > ^2 . 



(24) 



[ i.e. 6 = 0, in Eq. ], the length scaling of the gap crosses over to the stretched exponential form conjectured above 
in Eq. ([23]) . This is corroborated in the inset of Fig. [2] which for large sizes (and open boundaries ), supports 




FIG. 3: Non-universal dynamic exponents of the binary distribution of Fig.[2l Filled circles near pc stand for lower bounds of 
z. In the inset these latter support the divergence trend of the analytical results close to criticality [Eq. (|2ip in dashed curves]. 



the universal stretching factor (3 = 1/2 involved in (|23p . Thus, in nearing pc the stochastic dynamics slows down 
dramatically which, in line with the approach of the preceding section, is reflected in the abrupt increase of dynamic 
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exponents observed in Fig.[3l As we have seen in Eq. (|2T|) . for p ^ pc these latter can be arbitrarily large in the 
thermodynamic limit, but in practice they are ultimately cutoff by our available system sizes. In this regard, notice 
that the trend of slopes obtained near the critical region (e.g. lowermost curve of Fig. [2]), constitutes a numerical 
lower bound for these exponents (filled symbols of Fig. [3]). 
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FIG. 4: Distributions of gaps (g) for binary disorder over 500 sites. In (a) the parameter rates are taken as in Fig. [J] the 
main panel shows (part of) the case p — pc (o / {g) ~ 2), whereas the inset displays a characteristic non-critical situation {p — 
0.5, o-/{g) ~ 0.6). Both results have exponential tails (in bold lines), (b) shows a typical well defined bias case ({0.6,0.2}, p = 
0.5), exhibiting instead a narrow distribution closely following a Gaussian (in dashed line, o-/{g) ~ 0.02). As is shown in the 
inset, here dynamic exponents are difi^usive. 

Also, in agreement with the theoretical analysis, it is seen that the distribution of spectral gaps becomes broadest 
at Pc (maximum variance a). This is illustrated in Fig. 3^, after binning gaps in 800 intervals over 60.000 disorder 
realizations. By contrast, chains with non-fluctuating bias orientations have rather narrow gap distributions, and 
typically diffusive dynamic exponents (z — 2), regardless of the value of p. This characteristic scenario already 
anticipated in Sec. Ill, is exhibited in Fig.[4lD. 

Simulations.- To allow for an independent examination of nonuniversal aspects, we also carried out numerical 
simulations in periodic chains of 200, 400 and 800 sites (without duplication). Averages of particle densities were 
taken over 10'^ histories of 10"^ samples with binary disorder. Starting from homogeneous densities p, our results in 
Fig. El clearly indicate an asymptotic decay of the form (pL(t) ) oc e"*/"^^ with a relaxation time tl oc . As expected, 
the slope collapse follows closely the non-universal dynamic exponent already identified by direct diagonalization 
(z ~ 1.64, in Fig.[2])- In turn, the best scaling behavior for well oriented biases is obtained with standard diffusive 
exponents, which is also in agreement with our previous numerical results (inset of Fig.|3)3). 

Finally, it is interesting to comment on the regime i <C i^. This latter is addressed in Fig. [S] where the density decay 
below, at, and above Pc is exhibited for much larger systems. After averaging 20 samples over 100 histories each, it 
is observed that in all cases there is an early stage showing the typical decay. However, at p = Pc {p > Pc), an 

incipient crossover to a slower (faster) decay clearly emerges (cf. [E, ) • It would take further evolution decades to 
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<P(t)> 




0.5 1 

FIG. 5: Density decay for the distribution of Fig. [5] exhibiting a non-universal dynamic exponent at p = 0.5. The inset shows 
the case of a well defined sign of binary biases {0.6, 0.2} (as in Fig.|4l3), displaying instead a typical diffusive exponent. Squares 
denote data for L = 800, circles 400, and triangles 200. To compare slopes, the data of these latter two sizes were shifted 
downwards. 

estimate the actual form of such regimes, ahhough their nonuniversal characters are aheady in line with our findings 
for z(j)). The situation for p < pc remains unclear as the existence of a very slow crossover can not be ruled out. 




10 10^ 10^ 



steps 

FIG. 6: Density decay for the binary distribution considered above, using L — 10* and different concentrations. From top to 
bottom these vary as pc 0.2263 (circles), 0.15 (triangles) , and 0.5 (squares). For comparisons of large times trends, the data 
of the latter two concentrations were shifted downwards. Solid lines are fitted with slopes -1/2. 



V. CONCLUSIONS 

To summarise, after identifying phase steps, and limiting phases and their relationship to generalised biased random 
walks, in Sec. Ill we have outlined how the character (power law or stretched exponential) of the non-universal critical 
dynamics and in particular the critical exponents and their non- universality are determined by the associated diffusion 
constant and bias (in the case where the fundamental phase step can have either sign). 

Among the specific results given above are: in the general biased case the dynamic exponent z is -D/|6| and 
associated distributions of the energy gap are exponential. For vanishing bias z diverges and the usual critical power 
law behaviour is replaced by a stretched exponential form, possibly signaling the emergence of a glassy dynamic, 
while the ratio of rms gap to mean gap diverges. All those findings were corroborated numerically using the doubling 
fermion approach given in Sec. II and Monte Carlo simulations of Sec. IV. 

Concerning occurrence of both positive and negative phase steps, that is connected, in a given configuration of 
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randomness, to a succession of domains of alternating bias, whose boundary sites alternate between being traps or 
repellers for particles. That causes some particles to wander ineffectively for long times before finding a partner for 
pair annihilation (see Fig.[T]), leading to large or divergent dynamic exponents (see Figs. [2] and [3|. 

Among many important things not covered are: the full parameter-dependence of dynamic exponent z from pure 
value to divergence where the bias vanishes (e.g. in the case of the binary distribution of rates, with concentration 
parameter p, where z(p) is not monotonic for p in [0 , pc] or in \pc , 1]); the effect of open boundaries versus PBC (can 
be very significant in driven non-equilibrium systems) - particularly here regarding the stretching exponent for zero 
bias - needs further investigation. 
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